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Abstract. Doughnut shaped light beam has been generated from Gaussian 
mode (TEMqq) cw-Ti sapphire laser. After splitting the pump beam into two 
equal intensity components and introducing unequal convergence and phase delay 
while they are recombined it results in doughnut mode. Such a beam is tunable 
and have long propagation length. The evanescent field generated by 360 mW 
(at 780 nm wavelength) of such a beam creates optical field of 600 nm decay 
length with a 5.75 neV repulsive dipole potential. Thus cold Rb atoms (at 10/xK 
or less temperature) released from MOT can be reflected by the surface so that 
the atoms are collected ultimately at the bottom of the prism. By focussing such 
doughnut beam with 8 cm focal length converging lens, the dark radius reduces 
to 22fi. We also observe such beam to contain azimuthal phase as well as radial 
phase distribution. 



PACS numbers: 42.50.Vk, 42.60.Jf, 32.80.Ys 
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1. Introduction 

Localization of tiny particle by doughnut beam was first demonstrated by Ashkin 
[l]which initiated interest on TEMq\ laser beam and its application in atom 
manipulation. Besides, evanescent field atom mirror proposed by Cook and Hill [2] 
and subsequent interest on atom mirror led to observation of atom reflection [3] . In the 
same way confinement of cold atom by hollow core laser beam was realized [4] . These 
two mechanism of atom photon interaction have been nicely manifested through the 
proposal of atom funnel experiment [5]. Furthermore this hollow core beam (TEMqi) 
has drawn increasing attention because of its usefulness to manipulate and guide [6] 
cold atoms through its dark core, its ability to confine tiny particles [1] as well as atom 
reflection and collection in atom funnel experiment [5]. With computer generated 
hologram such a beam has been generated [7]. It is also possible to use cylindrical 
lens [8], or modified intracavity mirror of laser [9] to generate TEMoi mode. Using 
hollow core optical fiber a doughnut shaped optical near-field was also proposed [10] 
in order to guide cold atoms. In atom tunneling a double cone prism was used to 
generate cylindrical optical field distribution from a Gaussian beam [5]. The success 
of the Ashkin trap was the scattering force of photon generated by strong optical field 
distribution around the particle. However it gradually emerged that the theory of 
dipole force proposed by Askar'yan [11] can also help to control atoms [12]. For such 
application doughnut beam can be coupled to a hollow prism container receiving cold 
atoms so that localized surface waves are generated to collect cold atoms and preserve 
its state through reflection or, thereafter transport the atoms thought a bottom hole 
[5]. 

Various degree of dark-bright contrast and longer effective length of propagation 
are useful characteristics for its application in atom manipulation. With a presence 
of non-zero optical intensity at the center a longitudinal scattering force as well 
as transverse dipole force can effectively transport the atoms. Double cone prism 
produced beam [5] have steeper optical potential distribution across the bright ring, 
which is suitable for trapping, than that generated by hologram or other means, 
however such mode suffers quicker decay with propagation length. Here we report 
generation of various doughnut shaped beams that covers longer propagation length 
as well as various degrees of dark bright contrast and analyze its usefulness for atom 
manipulation. 

Cold Rb atoms are generated in a magneto-optical trap (MOT) . In order to collect 
a large number of such cold atoms a thin walled empty inverted pyramidal prism is 
kept below the MOT center so that after release from the trap the cold atoms fall 
into the prism under gravitational pull. When light have total internal reflection at 
glass-air interface an exponentially decaying evanescent field is created [2], a typical 
range normal to the surface is the optical wavelength of the light used. If frequency 
of the laser be slightly higher than the resonance absorption line of Rb atoms then 
the evanescent field can generate repulsive dipole force. Based on this principle blue 
detuned doughnut beam is applied through the bottom of the prism to induce enough 
atomic reflectivity at the prism wall and cold Rb atoms are collected by the prism [5] . 

2. Experimental 

In order to generate TEMqi beam we use two-beam interference in a Michaelson 
type interferometer. Figure 1(a) shows schematic diagram of the experimental setup. 
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A single mode Gaussian beam comes from Ti-saphire laser (maximum 1 W optical 
power). This beam is split into two equal intensity components by a non-polarizing 
beam splitter. These two beams are then passed through two different lenses F\, F 2 in 
order to focus along its path. Besides, to recombine and make these beams interfere, 
two mirrors are placed in a Michelson interferometer configuration so the beams are 
reflected back to the beam splitter. One of the mirrors have PZT actuator mounted 
at its back so that mirror displacement control in nm level be achievable. Thus the 
beams again pass through the respective lenses and have further convergence. The 
focal lengths of lenses F\ and F 2 are selected experimentally so that the beams are 
focussed close to the beam splitter after their second pass and their divergence beyond 
their respective focal points are appropriate to generate TEMqi mode. In this case 
the effective focal lengths are fx = 22.2 cm f 2 = 19.3 cm. This results in minimum 
beam waists of 110/i and 96/j, respectively at about 1.2 cm separation. Focal length 
of output collimator F3 is 40 cm. 




Figure 1. Doughnut beam generation (a)Schematic diagram of optical layout for 
doughnut beam generation. LIA indicates lock in amplifier, M is a small miror. 
(b) Topological diagram of beam convergence and collimation. 



2.1. Analysis of beam mode formation 

Each of the beams in the interferometer can be analyzed by taking Fourier transform 
of the input beam. This is because Fourier transform can convert frequency domain 
spectrum into spatial distribution, which means a monochromatic collimated laser 
beam while passed through a converging lens results in a single point in focal plane 
where it forms the image. 

Although doughnut beam has been described by Laguerre polynomial, Gaussian 
function, angular momentum quantum number [13,14], however it has also been 
observed that focussed Gaussian beam acquires a axial phase [15] which is known 
as Gouy phase. With this we try to analyze the result with simplified theoretical 
model specific to our system and look into various characteristics experimentally. 

Let E a (x,y,z) be initial electric field distribution (which is Gaussian in X — Y 
plane) of the optical beam (propagating along Z-direction), figure 1(b). For analytical 
simplicity let's start with one of the two beams. While passing through converging 
lens let the electric field amplitude be E(x' ,y' , z') at focal plane X' — Y' of lens F\, 
With the help of Kirchhoff integral: 



Doughnut beam. 



4 



E(x',y',z') 



E a (x, y, z)— exp[jfcr] cos(fc, f)dxdy (1) 

r 



r position vector, k wave vector, the angle (k, r) being very small cos(/c, f) ~ 1 , 
A optical wavelength (780 nm), j = V^T, initial beam waist as a . In this situation i 
can be expressed as ( Z /L z ) , and 
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With Fraunhofer approximation, exp [jjp- (x 2 + y 2 )] w 1 , however we take this 
factor outside the integral rather than replacing by unity. Thus by parameterizing 
constant factor in the above integral as A\ fi : 
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j jj- (a/ 2 + y' 2 ) is added due to the focal plane (X' — Y') , which 

if we compare to reference [15] has some relation to Gouy phase. Now we get from 
equation 1, for i = 1, 2 



E i (x',y',z')=A lfi 
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where fi focal lengths of lenses. With E a (x, y, z) = Aa exp — x ± , where wo 

L w o J 

is initial beam waist of input laser, Ao its electric field amplitude Ei (x',y',z') the 
integral can be written with the help of reference [16], as 



E i (x\y , ,z') = A lfi ^yfm$^exp 

This is Gaussian at (X' — Y') i planes along with a phase distribution function 
A\ f i . Its minimum beam waist as = Woi . However beyond the focal plane the 
beam diverges due to its geometric structure, so the beam waist will linearly increase 
along its propagation. Let be half angle of divergence then each beam waist can 
be expressed as Wi = (w 0i + z[ tanWj). 



(x> 2 + y' 2 ) 



(2) 
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It has been observed that the two beams interfere after their respective focal 
planes and the topology of generated pattern does not change while the resulting 
beam diverges. Thus, while a third lens (F 3 ) is used at some other distance apart in 
order to collimate the output beam with desired waist size it does not introduce any 
modification to beam pattern. It has also been observed that position of this third 
lens does not affect beam topology either. Let the lens F3 be placed at X" — Y" plane. 
The resulting beam diameter is 5 mm. 

As the beam widths are finite and unequal so the interference pattern will depend 
not only on relative axial delay between these two beams but also on its spatial optical 
field distribution. To simplify the analysis, let's introduce a phase part related to 
(5 mod A) as (sin 6*) factor to E 2 (x',y',z 2 ) and (cos#) factor to E\ (x',y',z'), where 
9=^(5 mod A) , then 

E (x", y", z") = E x (x\ y', z') cos 9 + E 2 (x', y', z' 2 ) sin 9 
or 



E(x",y",z") = ^^U /ie xp 

+ V7T— 4^1/2 ex P 



(x" 2 + y" 2 ) 



(woi + (z 2 +5) tan^i) 
' (x" 2 + y" 2 ) 



(u>oi + z 2 tan\E , 2 ) 



sm t 



(3) 



This is the simplified generating function for different shape of the output beam 
for delay 5. This is called 'simplified' because of a number of significant assumptions 
have been made so far. First and the prominent one is that the phase containing 
terms A\fi are equal for i=l and 2. Secondly we implicitly assumed that each of the 
two interfering beams do not have any azimuthal phase or helicity to alter intensity 
distribution of output beams. This is true for Gaussian beam however we will see 
later that while the output beam interfere with its mirror image, it exhibits transverse 
phase distribution. The possible beam shapes that are obtainable experimentally is 
shown in figure 2,3. 




(a) 



(b) 



(c) 



(d) 



Figure 2. Various doughnut beam images at output A for various delay, (a) 
Taking this as reference, (b) 100 nm delay, (c) 250 nm delay, (d) 390 nm delay. 



As seen in the figure 2&3 not all the possible configurations correspond to 
doughnut shape except for figures within the range of figure 2. As we see from the 
above expression (3) that, at x" = = y" or, at beam axis, the E {x" , y" , z") becomes 
zero only when #«37r/4or— 7r/4. However in this situation the optical intensity at 
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Figure 3. Beam profile for beams in A direction, (a)-(d) the delay correspond 
to that of figure 2(a)-(d), (e) 450 nm, (f) 550 nm, (g) 650 nm, (h) 750 nm delay. 



the off-axial points does not vanish identically. In order to get a general expression 
covering all off-axial point we expand the exponential terms and rearrange to get 
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mm. We can verify that in equation (4) while x" 
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9.21 x 10 4 r" 2 both in meter 2 where 
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all values of z" and are vanishingly small at the beam edge. Furthermore, the term 
containing exponential factor increases faster with r" than the other two terms, so it 
can further be approximated as 



E(x",y",z") = A' 1 



ht 2 



exp [- (ti - t 2 ) 



(5) 



which is the expression for radial amplitude distribution for first order LG beam 
(LGqi mode). Thus, with an odd multiple of n or (2n + l)7r axial phase difference 
between the two beams we get LGqi mode. Similarly, when the relative axial phase 
difference is 2mr maximum brightness at beam axis occurs with 



E b (x",y",z") = A[ 
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h + h 
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The figure 3(d) and 3(h) correspond to additional phase delay, introduced by the 
PZT actuator of figure 1(a). 

2.2. Topological Symmetry 

Thus we see from equation (3) that the achievable beam modes repeat after each A 
delay. The different beam profiles thus generated have topological symmetry, meaning 
the transformation from one beam structure to another is continuous, reversible and 
forms a closed set. Figure 3 shows the full range of available beam modes at several 
discrete optical delay in the range of 8 = 780 nm, and for any delay 6 more than 
integral multiple of A, the (6 mod A) parameter becomes effective and the generated 
beam mode falls within that shown in figure 3. Thus we can categorize the set (figure 
3) as a topological space with a pair (X, T) consisting of a set X of collection of possible 
beam modes and a family T of subsets of X so that: (a) an empty set $ e X and 
X e T (b) T is closed under arbitrary union (c) T is closed under finite intersection. 
Although topologically it is not allowed to pierce or join two ends of a structure because 
of possible ambiguity in regeneration and reversibility, however creating a dark or 
bright center, the set of beam profiles that can be regenerated are unambiguous and 
without any arbitration. This argument is valid for zero intensity points. In this we 
can introduce a blank screen as empty set (<J>) which while added or subtracted from 
any member does not introduce any change, in practice this corresponds to generating 
delay among the two beams so that a particular beam configuration can be maintained, 
which in other words means the beam mode locking. Embedding this theory into an 
ideal experimental situation it corresponds to zero delay to a stabilized configuration. 
The condition (b) states T is closed under arbitrary union, which implies that with 
any possible delay the achievable beam profile will always lie within that shown in 
the figure 3 or expressible by equation (3) as the effective delay (8 mod A) always lies 
within to 780 nm. Thus any subset of X are open and the collection T of open 
subsets can thus be called topology for X. 

With this categorical simplification we can further classify this experiment and its 
related results to look into its further characteristics. To proceed we broadly devidc 
the set X into two subsets, one with darker axial profile, figure 2 or figure 3(a)-3(d), 
and the other with brighter axial intensity distribution, figure 3(e)-3(h). Let these 
be subset A and B. As we know that in the Michaelson interferometer one can get 
two beams, one goes away from the laser and other return towards the laser. We 
identify this second beam with a prime e.g. B' for beams with brighter axial profile 
like that of figure 3(e) to figure 3(h) that returns towards the laser. We select A and 
B', as these are simultaneously observable complementary beam structures the sum of 
which always gives the total beam intensity (assuming zero loss through the optics). 
Similarly B and A' can form another complementary set. Thus the delay that results 
in evolution of A simultaneously generates B' . To look into this situation analytically 
we take the following logical and analytical method. 

A homeomorphism between spaces A and B' is a continuous map / : A — > B' 
which has an inverse continuous map g : B' —> A, so that f(g(b')) = b' for all b' £ B', 
and g(f(a)) = a for all a e A . Thus the maps / and g introduces one-to-one and onto 
function (bijection) between the spaces A and B' . We can see from experiment that A 
and B' are simultaneously generated and forms two subsets of topological spaces that 
are distinctly measurable and identifiable in real space. In the following we will look 
into explicit expression of / and g in accordance to above formalism and investigate 
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how this leads to interrelating the two spaces A and B' to stabilize A using information 
from B' . 

We formulate the / and g in terms of delay and represent with 9, which is a 
function of 5. For / and g as continuous and inverse of each other it is sufficient 
to show that the square matrix representing / and g are inverse of each other. Let 
cos 9 — sin 6 



/ = 



sin 9 cos 9 
and by Gauss Jordan method 

cos 9 sin 9 
- sin 9 cos 9 

thus g = Now rewriting expression (3) as 

E(x",y",z") = [pi cos9 + qi sin0] 

where 



and 



Pi = A\ exp 



qi = A\ exp 



V 



(w i + {z% + 5) tan* 2 ) z 



(x" 2 + y" 2 ) 



(w 02 + (4 + <5) tan fa)" 

represents the same physical situation in a convenient form. We can consider p x 
and q\ as bases of the topology because any member of T can be expressed as a linear 
combination of p\ and q\. Taking E [x" ,y" , z") = a and b' = E' (x" , y" , z") say, a 
matrix equation can be formed as follows 



a 
b' 



Pi 

qi 



cos 9 sin 9 
— sin 9 cos 9 



Pi 

qi 



(7) 



Or a = [pi cos 9 + q\ sin 9} and b' = [pi cos 9 — qi sin 9] . In practice the negative 
sign in b' arises as 9 should change by w for the beam B' . Thus the transformation 
relations corresponding to / and g are unitary, reversible and the elements that are 
generated in this way preserves total optical energy. Consequently, as is evident, 
for the delay (or 9) as independent variable of the transformation, the beam intensity 
corresponding to a and b' becomes unequal, except for 9 ^ (An + 1) 7r/4, which implies 
one of two beams acquires higher intensity in expense of the other. Table 1 shows the 
optical power of doughnut beam at various stages of beam mode. 

A point to be noted is that, although a and b' represents members of A and B' we 
can not use them in the above expressions to analyze their nature of continuity, that is 
test for monotonously increasing or decreasing nature will fail. This is because, as has 
been mentioned earlier, both a and b' are not uniform or their identity at different 9 
are not similar because of topological modifications. Thus, one can not get condition 
for maxima or minima by equating first derivative to zero. 



3. Results 

3.1. Doughnut beam modes 

It has been observed that the output power of the doughnut beam varies with mode 
of the output beam. Table 1 shows the variation of beam mode and their respective 
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Table 1. Optical power of doughnut beam A at various beam modes 



Laser power (mW) Beam mode Output (A) power a : b' 



250 Figure 2(a),3(a) 145 2.3:1 

Figure 2(c),3(c) 125 1.5:1 

Figure 2(d), 3(d) 105 1.0:1 



610 Figure 2(a),3(a) 360 2.3:1 

Figure 2(c),3(c) 300 1.4:1 

Figure 2(d),3(d) 260 1.0:1 



900 Figure 2 (a), 3 (a) 620 2.2:1 

Figure 2(c),3(c) 460 1.5:1 

Figure 2(d),3(d) 380 1.0:1 



optical power measured. Table 1 shows the optical beam intensity at various mode of 
the output beam [17, 18]. 

Note that there is some loss of light at the optics, however the measured A : B' 
beam power ratio remains constant. 

The figure 1 shows that the optical delay is controlled by piezo-electric actuator 
attached to one of Michaelson-mirrors, it controls the beam mode. Besides, there may 
be electrical noise, external mechanical vibration as well as some other disturbance 
may create unsteady situation unless and otherwise special precaution is adopted. We 
have already observed that the set of beam modes form a topological set where A and 
B' are interrelated (7) and the null element of which while added or subtracted from 
the set X does not have any effect. Now we shall use this idea to create a single element 
set or a stabilized doughnut beam output in which the null element will be functionally 
generated through locking circuit. This is useful in understanding locking the beam 
mode to a predetermined form. We know that / and g are inverse operations of each 
other. If there is some noise that arbitrarily generates / a complementary process g 
should be created to bring the system back to its initial state. We achieve this with 
the help of a photo-detector (PDi) and lock- in amplifier combination and feed this 
signal back to the PZT actuator (figure 1(a)) so that the feed back signal can stabilize 
the system. 

This has been achieved quite simply because the B' beam diameter is about 5 
mm. With a 0.5 mm at the PDi window aperture of PDi while placed in its path, can 
discriminate a minute change in beam mode. Chopping the beam B' and feeding the 
electrical output of the PDi to a lock in amplifier (LIA) the g is generated. This is 
more evident from figure 2 that beam intensity distribution continuously varies with 
radial distance from beam axis. Thus while a change of beam profile with rise in 9 
can lead to increased PD signal, a corresponding fall in 8 leads a change in PD signal 
on the reverse. 

3.2. Focussed doughnut beam 

One of several advantages that this beam have is its long propagation length; more 
than 3 meters observed in the laboratory. While doughnut beam is used to guide 
atom or as an instrument to manipulate atom in a localized space, the characteristics 
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Figure 4. Focussed doughnut beam (a)Flat-top aperture fiber probe as sampling 
window, (b) Measurement setup for minimum beam waist of free space focussed 
doughnut beam (/ = 8cm), (c) Doughnut beam profile measured at free space 
minimum beam waist by the probe. 
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Figure 5. Evanescent field, (a, b) protruded fiber probe used to detect optical 
near-field, (c)collimated doughnut beam in funnel-shaped evanescent setup, (d) 
evanescent field profile normal to the surface. 



of such focussed beam will be useful. To observe doughnut beam mode under focus 
we used gold coated flat-top aperture fiber probe (figure 4(a)) to pick up optical signal 
from minimum beam waist. 

Figure 4(b) shows the schematic diagram of free space measurement setup. The 
fiber probe is mounted on a 3-axis nanomover stage (Melles Griot Inc.) in combination 
with a PZT actuator of 60 nm/Volt sensitivity. Figure 4(c) shows the measured free 
space beam profile. Accuracy of free space measurement has some limitation due to 
scattering at the fiber probe as well as non-zero transmission of the gold coating that 
covers the probe tip. Single tapered fiber probe has been fabricated by conventional 
chemical etching process and making a 180 nm gold coating over the tip. A 370 
nm aperture has been created by focussed ion beam (FIB) technique, figure 4(a). A 
details of such method can be found in [19]. 

After that we carried out the optical near-field measurement where the optical 
field is "blocked" beyond the probe tip and the scattering effect becomes minimum. 
Figure 5(c) shows the schematic diagram of the optical near- field measurement setup. 
Figure 5(a,b) shows the probe used in this measurement. Figure 5(d) shows the optical 
near-field distribution of the focussed doughnut beam measured in this way. About 
20/i dark diameter has been recorded. It shows a dark hole probably exists along the 
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minimum beam waist at the focal plane and the TEM 01 mode can have a continuous 
dark core even though the focal plane, when it is focussed. Such a characteristic can 
be useful while transporting cold atoms through the dark core and manipulating in a 
small dimension. 

3.3. Evanescent optical near-field 

The focussed doughnut beam profile as discussed above, has been the result of optical 
nearfield measurement. So at this stage we will discuss about optical near-field 
generated by doughnut beam in atom funnel system, measure the evanescent field 
and estimate most probable behaviour of Rb atoms at the funnel wall. 

3.3.1. Evanescent atom funnel In order to measure the optical near-field generated at 
the funnel surface of the atom funnel experiment, we use bare long fiber probe, figure 
5(a), (b). A similar chemical etching technique, as described before, has been used 
with longer etching time to create long tip fiber probe for evanescent field measurement 
in funnel. Fiber probe-tip angle is about 23 degrees whereas for a short probe tip 
length of 6/i at fiber core of diameter 40/i. Thus in the evanescent field measurement 
over hollow prism it is inconvenient to use short fiber probe (for funnel surface, as the 
vertical wall angle is about 45°). As discussed earlier total internal reflection generates 
surface wave or evanescent field. The doughnut beam while is incident in such a way 
that a total internal reflection takes place, then a doughnut shaped evanescent field will 
be created at the prism surface. The measured optical near-field intensity distribution 
is as shown in figure 5(d). 

The fiber probe is placed within the evanescent field range and it is kept in a 
predetermined position by applying a small vibration to the probe tip and collecting 
corresponding error signal. 

Application of doughnut beam in atom funnel has been developing [5]. In this 
system the necessary beam profile should have at least 250/1 dark diameter along with 
2.5 mm beam waist. Cold atoms are generated in magneto-optical trap and collected 
by the hollow prism with the help of the evanescent field created by doughnut beam. 
The prism wall thickness is about 3 mm and is composed of three 2.5 x 2.5 cm 2 
glass plates glued together. The lower tip of the prism is polished to a flat triangular 
surface, at the center there is a 250/i hole. Doughnut beam is aligned at the bottom 
of the prism so that dark center lies at the prism hole and the beam gets total internal 
reflection at the inner surface. This illumination is synchronized with release of cold 
atoms from MOT so that the atoms bounce over the inner wall of the prism and 
collected through the bottom hole 

At three different optical power evanescent field of various strength have been 
measured. However the measured range of the field primarily depend on the detector 
noise level and collection efficiency of the fiber probe. Thus with intense doughnut 
beam we can detect longer range of the field. 

The characteristic decay length at which the evanescent field intensity decays to 
e _1 to that at the glass surface has been estimated from experimental result. The start 
of the exponential decay has been taken as the maximum value for the evanescent field. 
The minimum of the exponential decay corresponds to NEP of the detector system. 
Thus the range of the evanescent field is the range from the maximum evanescent field 
to the NEP level, which has also been estimated from experimental results. 
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The probe tip cross section has about 42 nm diameter. Here we can make another 
practical approximation. While the doughnut beam is incident from the bottom of 
prism, the three planar walls of the prism surfaces will be uniformly covered with 
optical near field . This is because the wall thickness is about 3.0 mm whereas each 
wall area is 2.5 x 2.5 cm 2 and the angle that the three edges make with vertical is 45 
degrees. Besides, width of the bright ring is about 2.5 mm. Combining all these with 
repeated total internal reflections both at the inner and outer surface of the walls we 
can approximate to an almost uniform distribution of optical near-field. Thus with 360 
mW beam the areal intensity distribution will be 213.3 W/m 2 . With 5.54xl0~ 15 m 2 
fiber probe area optical field will be collected and measured. Assuming 70% collection 
efficiency at the probe tip, because of refractrive index mismatch between probe and 
air, and also there be another 30% optical power loss at fiber-photodetector (PD 2 ) 
coupling. To calibrate optical field distribution inside the funnel wall we take the noise 
equivalent power (NEP) of the system as the reference. Observed LIA noise level is 
1.2/xV. With a 10 9 gain of PD 2 gain circuit, it corresponds to 4.1xl0~ 2 pA current 
generated at PD 2 . With a 50% photo- voltaic conversion, this corresponds to 2.4xl0~ 2 
pW optical power at the PD 2 , (although the NEP level of the PD 2 , at 200 Hz chopper 
frequency is 1.4xl0~ 4 pW). This correspond to 4.1xl0~ 2 pW evanescent light at 
the fiber probe tip. At 780 nm wavelength the photon energy is 1.59 eV, and thus 
the NEP corresponds to 9.4xl0 4 evanescent photons floating near the fiber probe tip. 
This is at a distance of 156 nm from the prism wall and at 360 mW doughnut beam 
power. Natural Rb atomic radius is 2.98A, so its cross sectional area is 2.8xl0~ 19 m 2 . 
Therefore, while Rb atom approach the prism wall normally the number of photons 
it come across is 5.0xl0~ 5 times the number of photons faced by the detector probe. 

3.3.2. Dipole force Grazing incidence reflection of atoms by evanescent field has 
been discussed since 1987 [20,21]. A blue detuned laser light while gets total internal 
reflection in glass gives a layer of repulsive evanescent field over the glass- air interface. 
At the blue detuning of laser line, because of faster change in laser field compared to 
that of atom it responds in a time frame lagging with respect to the external field. 
We can assume a frequency range (seen by the atom, or in atomic frame) in which 
this dipole force is appreciably high, lies close to resonance frequency of Rb atom. It 
is obvious that it can not happen at all higher frequencies or in any range of blue 
detuning. The principal difference between photon scattering (or photon reflection by 
atom) and this dipole force is that in the former case the incident photon frequency 
remains almost same to that of the scattered one, whereas in the latter its wave nature 
takes precedence. Thus without going into quantum mechanical analysis using dressed 
state formalism as well as drawing a sharp boundary between these two interaction 
mechanisms we can treat this atom kinetics semi-classically. 

In free space if we assume a stationary photon distribution that corresponds to 
total doughnut beam power then the Rb atom at 10/xK temperature will come across 
about 250 photons per second on its path, which in this case will be hypothetical 
limiting photon collision rate with cold and falling Rb atoms. To estimate the number 
of photons faced by Rb atoms at the prism wall, we take area under the " decay curve" . 
However one may argue that when Rb atom comes close to the glass surface a localized 
field (dielectric constant) will change and a large number of photons will tunnel across 
the barrier to collide with the Rb atoms. We eliminate this situation for ground 
state Rb atoms having low multipole moment, falling with small velocity normal to 
the surface, thus repulsive dipole force will be stronger before reaching the photon 
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tunneling range. Moreover, if Rb atoms do have sufficient velocity to reach to the 
photon tunneling range, the atoms will again be reflected due to velocity dissipation 
and photon scattering. Thus not only very slow atoms can reflect from the surface 
but also at higher atomic velocity the atom reflection can again be observed. 

The optical dipole potential on atom can be expressed in terms of detuning A as, 
Uo = TiAcr/2, the factor a k, (r /2A) 2 TI /Jo incorporates laser intensity I, saturation 
intensity of Rb atom Iq = 1.65 mWcm~ 2 , its natural linewidth T = 37.7 MHz, and 
a factor depending on polarization of the incident beam T« 5.7 (2.57) for TM(TE) 
polarization [22] and U = ?7o cxp [— 2 Act] where k = ^-n 2 \/sm 2 a — 1; n is refractive 
index of glass (with respect to air, 1.45) and a angle of incidence 45°. However, 
instead of using theoretical value for 2k, we shall estimate it from experimental result. 
Thus the retarding dipole potential is distributed functionally similar to optical field 
intensity distribution. With 0.36 W doughnut beam power, the dipole potential 
distribution can be written as Udip — 5.75 exp (— 2nt) neV, for A = 1.5GHz, where t 
is distance measured normal to the surface. 

Apart from atom optical interaction described above there will be attractive van 
der Waals force that acts on Rb atoms due to proximity of the glass surface, Uvdw(t) = 

~ 2 68(^+1) (kft) Li 1 [^1 where k = ^mT 1 . ft has been known that atoms at a distance 
of more than 0.2A the effective force is the repulsive dipole force and Van der Waals 
force have little effect on atom repulsion beyond such a distance [12]. Thus considering 
the reflected atoms we can neglect the Van der Waals force at a distance longer than 
156 nm. So the repulsive optical force can be written as 

U op t ~ 5.75 exp (— 2/tt) neV, for t >156 nm. We know from evanescent field 
measurement that 1/2/t = 600 nm. By integrating this potential from infinity to 156 
nm one can get work done on Rb atom while it moves toward the evanescent surface. 
Cold Rb atoms at 10/iK temperature have most probable kinetic energy of 1.3 neV, 
which is smaller than the 5.75 neV. Thus cold Rb atoms will be reflected by the 
evanescent field created by the doughnut beam at the funnel surface. 



3.3.3. Dissipation In addition to dipole kinetics near the evanescent surface, atom 
may also have dissipation of its energy at the evanescent field interaction [23] . Now we 
will consider dissipative effect of evanescent field interaction in which atom looses part 
of its velocity near the surface. Rb atoms (say 85 Rb) have two ground levels F=2,3 
separated by 3.036 GHz frequency. It has been observed that while blue detuning 
increases from 1 GHz to 1.5 GHz the number of reflected atoms also increases [24]. 
With this trend we can assume that due to atomic velocity and the consequent Doppler 
shift a transition takes place from F=2 state by absorbing the photon (frequency v\) 
and a subsequent emission at frequency v 2 to bring the atom back to F=2 state. Thus 
as the blue detuning increases so the transition probability rises. So from conservation 
of momentum, mu 2 — mu\ = (yi — v 2 ) h/c, where u\_ 2 are initial and final velocities 
respectively, u\ and v 2 are in the same direction. And from conservation of energy, 
Ei + hv\ + 0.5muf = E 2 + hv 2 + 0.5mu 2 . Combining these two w 2 — = m \t ran , ' 
where m is mass of Rb atom, and Xtrans corresponds to blue detuning (with respect 
to 5S F — 3 state), which in turn implies higher the blue detuning higher will be 
the atom reflection from the prism wall, as has been observed. Additionally we are 
considering atoms in motion so the atom will see the laser line with Doppler shift, 
which is AXoop — ±(wiA/c), we shall take negative sign while atom is moving towards 
light. So u 2 — mi = — > 7^ rr 7 which is the magnitude of velocity rise opposite 
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to the direction of U\. Numerically we get u 2 — U\ = 23nm/sec. It has been observed 
that evanescent field may act dissipatively to incoming cold Rb atoms [23]. 



3.3.4- Radiation force In addition to the dipole interaction, van der Waals force and 
velocity dissipation, there will be photon scattering from atoms as the atoms follow 
Maxwell Boltzmann velocity distribution. Let Rb atomic radius be r = 0.298nm. 
Then most probable number of atomic collision with an evanescent photon will 
be n — nr 2 l f ; 156 sq exp (— 2/ct) dt, where I is range of evanescent field, s =4.7; 
photon calibration constant. This can be stated in terms of the figure 5(d) as 
n = ittqIsoX (area under the plot). Numerically n = 1575 photons. The mechanism 
of scattering or photon reflection by atom can be described simply as follows. 
Conservation of momentum gives m(u2 — Ui) = 2/i/A or U2 — u\ = 12mm/sec for 
each head-on colision with one photon. The estimated temperature of Rb atoms 
falling on the prism is 10/iK [5]. Accordingly, the most probable velocity of Rb atoms 
will be 44 mm/sec. Thus the atoms need 4 'pure' scattering with photon if there be 
no induced dipole force and no velocity dissipation. As we can see from experiment 
the probable number of atom-photon collision is more than 6 (rather it is estimated 
to be around 1575), so atom reflection probability is high at 360 mW doughnut beam 
power. 



4. Azimuthal phase 

A characteristic property doughnut beam is its helicity or axial phase distribution 
[13]. Moving along circumference of such a beam cross section one faces varying optical 
phase covered by the path (or line integral of azimuthal phase) , which can be expressed 
with a factor exp(im<p) in its amplitude expression [7,13,14], where m is azimuthal 
quantum number. Moreover, as we can see from Kirchhoff integral that the factor ^4i/j 
lies outside the integral and we keep it as constant factor for integration. However, 
considering exp j^- {x 2 + y 2 ) factor of A\fi we can say that with r' — y/ x' 2 + y' 2 

and r = y/x 2 +y 2 of the beam its phase has additional parabolic distribution with its 
Gouy contribution [15]. Thus a radial phase (f) r = ^y- (r' 2 + r 2 ) is not constant over 
minimum beam waist and its nature of variation remains unchanged while it diverges 
beyond the focal plane. 

So the doughnut beam acquires both azimuthal [13] as well as radial phase 
distribution. To observe such a characteristic nature of the wavefront we performed 
interference of the doughnut beam with its self image. Figure 6(a) shows schematic 
diagram of such experimental setup. A 2.5 mm thick parallel face plane glass slide has 
been chosen for creating such a self image and interference by division of amplitude. 
As we can not modify separation of the faces of the glass slide so we adopt to change 
incidence angle (by several milli-radian) of the beam by rotating the slide in a vertical 
axis. At this stage we are not going to have quantitative analysis, so we will not 
discuss precise value of this angle and corresponding change in interference pattern. 
Thus a qualitative estimation will be sufficient for the discussion that follows. 

It may be noted that the radial phase distribution that arises within Kirchhoff 
integral is inherent to each of the two beams. Thus it is independent of this doughnut 
beam generation and will remain present in all beam modes of figure 3. Now we try 
to establish the fact experimentally that such phase distribution indeed is present in 



Doughnut beam. 



15 




(a) < e > (f) ® 

Figure 6. Doughnut beam interference with its self image, (a) Optical setup, (b)- 
(f) various dark-bright intensity distribution observed by CCD camera, (g) This 
figure demonstrates presence of helical wave front while other phase contribution 
is minimum. 

the resulting beams described above. Both transmitted beam and reflected beams of 
figure 6(a) have been recorded simultaneously, as shown in figure 6(b) - 6(f). 

In figure 6(b)-6(f) we see that interference dark lines, arising due to parallel 
surface slide glass is linear in some part while non-linear at some other. Linear dark 
fringes arise when optical phase relationship between the two beams is linear, and 
similarly non-linear dark line correspond to non-linear phase relation. Due to helical 
wavefront and radial phase distribution, it is possible in the experiment of fig 6(a) 
that phase of a part of wavefront becomes linearly related while the other turns out to 
be non-linearly related. Thus, while incidence angle is changed (in milli-radian range) 
from 6(b) to 6(f), a transition from linear to non-linear dark fringe takes place. 

The transmitted part of the beam, going through the glass slide, does not 
show repetition of such pattern because of intensity mismatch of component beams. 
Furthermore, the azimuthal phase distribution is also demonstrated in figure 6(g), 
which shows only one 'spoke' [13] of the interference image that corresponds to helicity 
m = \. 

5. Discussion 

We see that the approximated analytical expression of the resultant beam of 
interference is followed by the experiments performed, however the effect of radial 
phase and azimuthal phase distribution that we did not take into account at doughnut 
beam formation stage, and it becomes evident through later part of experiment. In 
figure 2, 3 we see that the dark core of the doughnut beam does not arise only when 
condition 9 = (An+ l)w/4 is fulfilled, rather it tends to prevail through a range of 
delay. It may be due to radial as well as azimuthal and Gouy phase distribution of the 
two interfering beams dominate in such a zone which results in slow transformation 
of the dark axial region to bright zone. 

In case of evanescent field and its role on atom reflection we looked into three 
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possible mechanisms of atom photon interaction among these the dipole force and 
retardation effect of the optical near-field are polarization dependent while the 
scattering force does not depend on beam polarization. In addition dipole force is 
weaker than scattering force. Therefore, if we extend the analysis a bit further, it may 
be possible that the scattering force and the retardation effect have greater role to 
play in atom guidance of cold atoms through dark core laser beam. 

Van der Waals force is the force in atomic range effective among atoms. When Rb 
atom comes closer to the glass wall its electric multipole moment and that of the glass 
surface will be responsible to generate such attractive force. The electric multipole 
moment will depend on atomic state, e.g. in the ground state the moment will be 
lower than that at excited state. However atom in the excited state while enter the 
evanescent region may subject to emit a photon in presence of dipole force field and 
get de-excited. As a result the atom may face a varying strength of Van der Walls 
force. However under the condition when Van der Waals force is weaker than dipole 
force and dipole force weaker than scattering force, the effect of such a localized and 
instantaneous effects may not have major effect in atom reflection in the atom funnel. 

To transport cold atoms through its dark core it is necessary to have a suitable 
handle so that preferential force on such atom can be applied. A very weak but non- 
zero axial optical intensity profile can provide such a force through scattering. This in 
combination with radial optical potential field of the doughnut beam may transport 
and manipulate cold atoms more efficiently. 

6. Conclusions 

Thus, we have generated tunable and stabilized doughnut beam through interference 
in a modified Michael interferometer. The optical power is sufficient to reflect cold 
Rb atoms, at around 10^K temperature. It has been observed that such a beam have 
radial as well as azimuthal phase distribution and helicity as i . Out of various modes 
that available in this method, a doughnut mode have beam waist of 2.5 mm and dark 
diameter more than 250fi. 
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